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O ■ Abstract 
(N ■ 

' An explicitly energy-conserving full orbit code CUEBIT, developed originally 

D ■ to describe energetic particle effects in laboratory fusion experiments, has been 

applied to the problem of proton acceleration in solar flares. The model fields 
^—^ I are obtained from solutions of the linearised MHD equations for reconnecting 

■ modes at an X-type neutral point, with the additional ingredient of a longitudinal 

^ . magnetic field component. To accelerate protons to the highest observed energies 

^ I on flare timescales, it is necessary to invoke anomalous resistivity in the MHD 

solution. It is shown that the addition of a longitudinal field component greatly 
ff^ . increases the efficiency of ion acceleration, essentially because it greatly reduces 

the magnitude of drift motions away from the vicinity of the X-point, where the 
accelerating component of the electric field is largest. Using plasma parameters 
• consistent with flare observations, we obtain proton distributions extending up to 

I 7-r ay-emitting energies (> 1 MeV). In some cases the energy distributions exhibit 

Qh' a bump-on-tail in the MeV range. In general, the shape of the distribution is 

Q ■ sensitive to the model parameters. 

^: 
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The process of magnetic reconnection is believed to be intrinsic to solar flares. It is 
invoked as the mechanism whereby energy-loaded magnetic fields can reconfigure to 
a lower energy state, liberating the energy which powers particle acceleration, heat- 
ing, and mass motions. As the resistivity of solar coronal material is very low, the 
large reconnection rate required to power a flare demands the presence of reconnect- 
ing structures with very small length scales - namely, current sheets. Both in solar 
and magnetospheric physics, the study of particle acceleration in such structures has 
received considerable attention over the past few decades. Under non-steady state 
conditions, a reconnecting magnetic field produces an inductive electric field E: the re- 
sistive magnetohydrodynamic (MHD) form of Ohm's law indicates that in general this 
E has a component parallel to the local magnetic field, and can thus accelerate charged 
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particles. An example of such a non-steady state is a perturbed reconnecting X-type 
neutral point, for which Craig and McClymont (1991) calculated the normal modes in 
a strictly 2-D field, B{x,y). In this paper we demonstrate that the analysis of Craig 
and McClymont is applicable when a ^-invariant component of the magnetic field is 
included, enabling it to be extended to configurations more likely to be representative 
of conditions in the flaring solar corona. 

In the context of solar flare physics, a current sheet is often envisaged as arising 
in a two-ribbon flare, where reconncction of oppositely-directed, predominantly verti- 
cal magnetic fields results in a sheet structure having a length (parallel to the solar 
surface) on the order of the length of a post-flare arcade (~ lO^m), a vertical extent 
(height perpendicular to the solar surface) comparable to this, and a thickness on the 
order of a few times the ion gyroradius. Particle acceleration - in particular proton 
acceleration - in just such a (collisionless) reconnecting current sheet was addressed by 
Martens (1998). Following from the magnetotail work of Speiser (1965) he considered 
the effects of the inclusion of a component of magnetic fleld perpendicular to the plane 
of the current sheet, deriving expressions for the maximum energy attained before the 
gyromotion around this perpendicular fleld causes the protons to exit the sheet, and 
hence the acceleration region. With reconncction electric field strengths comparable 
with those observed (e.g. Kopp and Poletto, 1986), it was found that the proton energy 
and flux budgets for flares were compatible with those arising from acceleration in a 
macroscopic current sheet (for example formed in the wake of a rising fllament, with 
vertical extent on the order of lO^m, comparable to its length). 

Litvinenko and Somov (1993) also considered the current sheet geometry, demon- 
strating the effects of the addition of a further magnetic field component parallel to 
the direction of the reconncction electric field. A component of the charged particles' 
motion is then gyration about this parallel field, which allows them to stay longer in the 
current sheet, and attain higher energies. It is this fact which also becomes important 
in our studies. 

While a reconnecting current sheet with finite (and large) vertical extent is fre- 
quently used in considering acceleration in large two-ribbon fiarcs, where such a con- 
figuration is expected on the basis of numerical simulations of filament lift-off, we 
concern ourselves here with a current sheet of zero vertical extent, i.e. a coronal X- 
line, which corresponds to the initial stages of reconncction behind a filament. In the 
simplest case, this configuration is identical in cross-section [the {x, y) plane] to a 2-D 
X-type neutral point, but is invariant in the ^-direction. In reality, true magnetic nulls 
[B — 0) are likely to be rare in the corona. The elaboration of this scenario which we 
study here, and find to have important consequences for particle acceleration, is the 
inclusion of a component of magnetic fleld also in the ^-direction. 

Proton acceleration at an X-point with flnite was investigated by Bulanov 
(1980), Mori, Sakai, and Zhao (1998), Browning and Vekstein (2001), and Bruhwiler 
and Zweibel (1992). The present study differs from that of previous workers in that the 
accelerating electric fleld Ez is obtained self-consistently from a magnetic flux function 
corresponding to a reconnecting eigenmode of the X-point. This perturbation solution 
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was invoked in the context of particle acceleration also by Petkaki and MacKinnon 
(1997), but without the addition of the finite component. 

A major challenge for any flare acceleration mechanism is to obtain sufficient high 
energy particles in a short timcscale; new observations from the RHESSI satellite are 
putting ever tighter constraints on timescales, spectra and total energies. The typical 
requirements for protons, as summarised by Miller (1998) are as follows: they are 
accelerated up to energies of ~ 100 MeV on timescales of about a second, and to about 
a GeV on timescales of a few seconds. Proton acceleration lasts for several tens of 
seconds, at a rate of ~ lO^^s"^, such that the total energy content in protons above an 
MeV is ~ 10^^ J. Given a coronal density of 10^^ — lO^^m"^ and volume of ~ lO^^m"^, 
this implies that each second around 1-10% of all coronal protons must be accelerated 
to MeV energies (and therefore require rapid replenishing - we do not address this 
here). We find that the inclusion of a moderate longitudinal field component greatly 
assists in this process. 

This paper is structured in the following way. In Section 2 we describe our new 
algorithm for particle calculations. Section 3 describes the model employed in the 
simulations. Some sample simulation results and the effect of varying parameters of 
the simulation, such as the strength of the z component of the magnetic field, are 
described in Section 4 and we end with discussions and conclusions in Section 5. 



2 Energy— Conserving Algorithm 

The nonrelativistic Lorentz force equations 

m— = Zev X B(x) + ZeE(x), -i- = v, (1) 

(Jjv (JjL 

are approximated in the CUEBIT (CUlham Energy-conserving orBIT) code by the 
following finite difference equations 

m =Ze{ xB ) + ZeE i , (2) 



At 



•jj-j+l -j^i Y*+l -|- 



At 2 • 

Here m, Ze denote particle mass and charge. At the start of each timestep x*+^ is set 
equal to x* in a first calculation of v*"*"^. A few iterations are made which converge 
quickly to a final v*"*"-*^. 

Mori, Sakai and Zhao (1998) also used a scheme in which v was set equal to 
(v*+^ _l_ v*)/2 on the right hand side of the Lortentz force equation. In the special 
case where E = 0, the scalar product of the right hand side of Equation (2) with 
v*"*"^ _l_ is identically zero, so that (v^^^Y = (f*)^. The scheme thus conserves en- 
ergy exactly. This makes it possible to obtain accurate results with relatively long 
timesteps. A modified version of Equation (2) conserves total energy exactly when 
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E is a finite potential electric field. For non-potential electric fields, such as those 
arising from magnetic reconnection, we have found the method remains accurate for 
large timesteps (exceeding the Larmor period in most of the computational domain), 
allowing large numbers of particles to be simulated for timescales that are relevant 
for flare acceleration. CUEBIT has also been benchmarked by using it to compute 
energetic particle orbits in magnetic fusion experiments (Wilson et al., 2002), and will 
in the future be used to study particle transport under various conditions (e.g. in the 
presence of turbulent electromagnetic fields) in such experiments. 

For the simulations discussed in Section 4 it is not necessary to incorporate rela- 
tivistic kinematics in Equation (2) as the energies reached are only a small fraction of 
the particles' rest mass energy (we consider only proton acceleration). However, a rel- 
ativistic version of the code is currently being developed for the purpose of describing 
electron acceleration and transport. 



3 Model of a reconnecting X-type structure 
3.1 Craig and McClymont solution 

For a simple model of a reconnecting field at an X-type neutral point we use the two- 
dimensional description of Craig and McClymont (1991), with the additional element 
of a finite third magnetic field component B^. To determine the conditions under 
which the Craig and McClymont analysis applies with Bz 7^ 0, we re-derive below 
their equations for the evolution of a fiux function ip, defined such that the curl of 
A = is equal to the magnetic field in the {x, y) plane. The induction equation 

^ = V X (v X B) + ^V^B, (4) 

can be written in the form 

^ = V X (V X A) + ^V'A. (5) 

ot HQ 

Here v is fiow velocity, 77 is resistivity (assumed to be constant) and jiQ is the per- 
meability of free space. The first term on the right hand side of Equation (5) can be 
expanded to give 

V X (V X A) = V (v ■ A) - (v • V) A - A X (V X v) - (A • V) V. (6) 

We assume that flows only occur in the [x. y) plane and that there are no variations in 
the ^-direction. In these circumstances the flrst, third and fourth terms on the right 
hand side of Equation (6) are all zero and Equation (5) reduces to 

^ + (v-V)V' = ^V^V', (7) 

01 Hq 
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which is the form of the induction equation used by Craig and McClymont. 

Neglecting plasma pressure, viscosity and any external forces such as gravity the 
momentum equation is 

|^ + (vV)v= Vb, (8) 

where j, p denote current and mass density. Now in the present case 

j - —V X B = -— VV- (9) 

Substituting this expression into Equation (8) and noting that V^' is in the (x, y) plane, 
we obtain 

+ (v • V)v = (VV) VV', (10) 

at ij,qp ^ ^ 

i.e. Equation (2.3) in Craig and MyClymont (1991). 



3.2 Extension of Craig and McClymont solution to 5^ 7^ 

We now assume that the B field is of the form 

B = V X (ijjz) + B,z, (11) 

where is constant and uniform. As before, we assume that there are no variations 
in the z direction. Substituting Equation (11) into Equation (4) and writing A — ipz 
as before gives 

d 

— (V X A + B,z) = Vx(vx(VxA + B,z)) + — (V x A + B,z) . (12) 
ot po 

Since B^ is constant and uniform this reduces to 

— (V X A) = V X (v X (V X A)) + V X [v X {B,z)] + (V x A) . (13) 

ot po 

This is equivalent to the induction equation of Craig and McClymont except for the 
term containing B^ on the right hand side. This term can be written in the form 

V X [v X (S,z)] = - (S,z) V ■ V. (14) 

Since B^ is not necessarily zero, the only way for this quantity to be zero is if V • v = 0, 
i.e. the plasma flow must be incompressible. This assumption is often invoked in 
studies of magnetic reconnection (e.g. Biskamp, 2000; Priest and Forbes, 2001), and 
is a reasonable one for solar flare plasmas. Incompressibility is implicit in the model 
of Craig and McClymont: they use Equation (10) in a dimensionless form that is only 
valid if the plasma density is constrained to be uniform and time-independent. In the 
limit of ideal MHD, it is straightforward to show that the flow must be incompressible 
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when, as in the scenario considered here, there are no variations in the 2;-direction and 
Bz is finite (Strauss, 1976). With the assumption of incompressibihty. Equation (7) is 
vahd for finite B^. Moreover, since Bz is curl-free it does not contribute to the current; 
since the latter is oriented in the z direction [Equation (9)], it follows that Bz does 
not contribute to the Lorcntz force either, and thus the momentum equation of Craig 
and McClymont [Equation (10)] also remains vahd. The solutions of these equations 
derived by Craig and McClymont are therefore applicable when the magnetic field has 
a uniform and static z component. 

3.3 Perturbation solution 

The equilibrium X-type neutral point structure invoked by Craig and McClymont 
(1991) is described by 

B = |^(yx + xy), (15) 

where Bg is the field strength at the (circular) boundary of the system where R = Rg, 
R being the radial co-ordinate measured from the z-axis outwards. Perturbations to 
this field geometry will result in normal modes of oscillation. Following Petkaki and 
MacKinnon (1997), we consider only the fundamental mode, with zero azimuthal and 
radial mode numbers. Approximate analytical solutions for the perturbed magnetic 
and electric fields can be obtained by dividing the system into an ideal outer region 
and a resistive inner region, separated by a critical radius 

i?e = i?o(|)\ (16) 

where S, the Lundquist number at the boundary R = Rg, is equal to iiqRoVa/t] where 
va is the Alfven speed &X, R = Rq . The solution for the perturbed magnetic field in 
cylindrical polar coordinates can be approximated by (Craig and McClymont, 1991) 

6B — —aB—ucos yum — j cos \^~^~^ j j ^ ' 

for R > Re- Here dB is an arbitrary scaling parameter, o; is a dimensionless mode 
frequency, a is a dimensionless decay constant and (f is the unit vector in the azimuthal 
direction. For R < R^ the magnetic field perturbation is zero. The value of uj is 
determined by the requirement that Equation (17) yields 5B = at i? = Re-, and 
a ~ ci;^/2 (Craig and McClymont, 1991). The electric field resulting from this solution 
is wholly in the 2;-direction and is given by 

Ez — dBvA sin In cxp (^~~^^^ (^^ ^ ^ ("^~^)) ' ^^^^ 

for R> Rc and 

7-, ,7-. / Ct'^A \ / f^VA \ f^VA W , 

= -dBvA^W y-^t j [a COS j +o^sm jj , (19) 
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for R < Re- Equations (17), (18), and (19) remain valid for an incompressible plasma 
with finite constant in time and space. We show below that the inclusion of a modest 
Bz dramatically increases the efficiency with which test particles are accelerated in this 
magnetic geometry. 

4 Acceleration in a reconnecting X-type structure 

4.1 Choice of parameters 

We consider specifically the acceleration of test particle protons in a prescribed field 
of the type discussed in the previous section. It is known that protons are accelerated 
in flares to energies of at least several tens of MeV in a timescale of the order of 
one second (Miller, 1998; Aschwanden, 2002). A straightforward calculation indicates 
that this could be achieved with a parallel electric field E\\ of around 1 Vm~^. Such 
field strengths are, incidentally, consistent with the parallel electric fields implied by 
observations of the separation rates of fiare ribbons, e.g. Kopp and Poletto (1986). 
However, as noted recently by Craig and Litvinenko (2002), fields of this magnitude 
are not consistent with classical Spitzer resistivity. In deriving their reconnecting field 
solutions, Craig and McClymont (1991) used the resistive MHD form of Ohm's law, i.e. 

= ?7J[| where jy is the parallel component of j. The maximum possible current density 
in a hydrogen plasma is jmax = 2nec, where n is particle density and c is the speed 
of light. So a coronal density of 10^^ — 10^^ m~^ implies that jmax ~ 10^ — 10^ Am~^. 
However, if the resistivity were determined by classical electron-ion collisions, and the 
plasma temperature T were equal to 10^ K, a field of 1 Vm"^ would produce a steady- 
state current density of around 2 x lO^Am"^. This high current cannot be sustained 
by a classical collisional resistivity. We deduce that protons can only be accelerated 
to 7-ray-emitting energies in a reconnecting field of the type invoked by Craig and 
McClymont if rj is much higher than the Spitzer value, i.e. it is anomalous. Craig and 
Litvinenko (2002) have pointed out another reason for ruling out Spitzer resistivity in 
the fiaring corona, namely that it implies a resistive scale length that is smaller than 
the collisional mean free path. 

Anomalous resistivity can result from lower hybrid or ion acoustic micro-turbulence. 
Whether lower hybrid or ion acoustic waves are more likely to produce the required 
anomalous resistivity depends on several parameters, including the ratio of plasma pres- 
sure to magnetic field pressure, the electron to ion temperature ratio, and the direction 
of the current relative to the local magnetic field (Aparicio et al., 1998; Biskamp, 2000). 
In either case, turbulence results from the current density exceeding a certain threshold 
KneCs where Cg is the sound speed and k is a numerical factor typically ranging from 
unity up to about the square root of the ion to electron mass ratio (Kulsrud, 1998; 
Aparicio, et al. 1998). The anomalous resistivity then prevents the current density 
from exceeding this threshold: the turbulence increases the effective collisionality of 
the plasma. In the case of ion acoustic turbulence, for example, there is a strong in- 
teraction with protons whose velocity component v parallel to the wave propagation 
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direction is of the order of Cg (Ishihara and Hirose, 1981). The strength of the interac- 
tion falls off roughly as and is neghgibly small for particles lying sufficiently far 
out in the tail of a Maxwellian distribution (cf. Litvinenko and Somov, 1993). There- 
fore, some small fraction of the initial proton population is effectively collisionless, and 
can be described using Equation (1): it is this sub-population of protons which we 
simulate using CUEBIT. The bulk of the plasma is effectively coUisional, due to the 
anomalous resistivity. The use of a test particle approach requires that any net current 
associated with accelerated particles is small compared to the current corresponding 
to the reconnecting field. 

For simplicity, the particles in most of our simulations were given zero initial veloc- 
ity. This is permissible, since the thermal spread of velocities corresponding to a typical 
coronal temperature is very small compared to the near-relativistic speeds required for 
protons to excite 7-ray emission, and we will show that that the final proton energy 
spectrum above ~ lOkeV in a typical simulation does not change significantly when 
the 5-function initial velocity distribution is replaced with a lO^K Maxwellian. 

The choice of parameters in the simulations is dictated in part by the hmits imposed 
by anomalous resistivity, combined with the typical parallel field strengths needed. We 
assume, for dcfinitcncss n = lO^^m"'^, Ro = lO^m and B„ = 0.01 T. If the temperature 
is assumed to be lO^K, and the current is assumed to be limited to necs (i.e. k, = 1), the 
anomalous resistivity required for E\\ = 1 Vm~^ is then 1.4x 10~^ flm and the Lundquist 
number at the boundary R = Ro {a. key parameter in the Craig and McClymont model) 
is ,5 ~ 2 X 10^°. Equations (3.6) and (3.8) of Craig and McClymont (1991) then yield 
a; ~ 0.13 and a ~ 0.0088. The scaling parameter dB in Equations (17), (18) and (19) 
was chosen to give Ez — 1 Vm~^ at t = and R < Rc (in this region Ez — E\\) rising 
to approximately 1.4Vm~l after Is. For comparison, simulations were also carried 
out with Ez = 0.1 Vm~^ and 10 Vm^^. Strictly speaking, a change in Ez implies a 
change in 77 and hence S if the limiting current is assumed to be a fixed multiple of 
neCg. However, since cu and a have only a weak (logarithmic) dependence on S, the 
same values of these parameters were used in all the simulations. 

The only remaining parameter to fix is Bz- It follows from the conclusions drawn 
in Section 3.2 that this can be varied freely without any of the other parameters being 
affected. Simulations were carried out with = OT, 10"^ T, 10"^ T, IQ-^^T and 
10~^ T: these are all reasonable values for active region coronal magnetic fields. 

4.2 Results 

In each simulation the trajectories were computed of approximately 10^ protons with an 
initially uniform random distribution of positions in the reconnection region {R < R^,), 
and, unless otherwise stated, a 5-function velocity distribution. The duration of each 
simulation, 1 s, was chosen on the basis that protons appear to be accelerated to tens 
of MeV on timescales of this order (Miller, 1998; Aschwanden, 2002). The timestep 
At was of the order of 100 Larmor periods calculated R — Ro (simulations with 
shorter timesteps produced essentially identical results); there were approximately 1500 
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timesteps per simulation. The value of dB was chosen such that the magnitude of the 
magnetic field perturbation a,t R — Rg and t — was always < 10~^ T. 

Figure 1 shows R and particle energy E as functions of time for a proton initially 
lying close to the X-point {R = 10m) in field configurations with Bz = (solid curves) 
and Bz = 10~'^T (dashed curves). The other field parameters are those listed at the 
end of the previous section. The proton initially lies at azimuthal angle ip — 24° 
in the {x, y) plane, and has velocity components in the cylindrical coordinate system 
Vr — — Vz = 10^ ms~^. 

It is immediately clear that the addition of even a very modest Bz has a dramatic 
effect on the trajectory of the particle in phase space. When Bz = the total magnetic 
field is very small close to the X-point (since SB — Oin this region). The Lorentz force in 
the {x, y) plane is consequently very weak, the particle is effectively unmagnetised, and 
hence moves rapidly away from the X-point. There is only a short phase of acceleration 
(in this case to about 1 keV): far from the X-point, the particle becomes magnetized 
and, since B is strictly perpendicular to E, there is no further acceleration. However, 
a longitudinal field of 10~^T is sufficient to confine the particle to within a few tens 
of metres of the X-point, with the result that acceleration to energies of 100 keV and 
above can easily occur. Because B^ is assumed to be uniform and the field components 
in the {x, y) plane are very small, the combined effect of grad-S, curvature and E x B 
drifts on the particle trajectory is negligible. 




0.0 0.1 0.2 0.5 O.'i 0.5 0.6 10 20 30 AO 

t (ms) t (ms) 



Fig. 1. Distance from the X-point versus time for particles initially 
lying close to i? = in field configurations with B^ — (sohd curve) 
and Bz = 10~^ T (dashed curve), (b) Energy versus time over a longer 
timescale for the particles whose drifts are shown in (a). 

Figure 2 shows final proton energy spectra for B^ = 10~''T and initial electric field 
in the inner resistive region Eqz — lVm~^. Two different initial proton distributions 
were used: a 5-function (faint curve) and a 10 ''K Maxwellian (bold curve). In both 
cases, protons are accelerated to energies of up to several MeV. The only significant 
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differences between the two distributions occur at energies < lOkeV: above this energy, 
the final distribution is insensitive to the initial conditions. 
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Fig. 2. Final proton energy distributions for Bz = 10 ^T, Eo^ — 
lVm~^. The initial distributions have temperature T = (faint curve) 
and T = lO^K (bold curve). 

Figures 3 and 4 show final proton energy distributions for several pairs of values 
of Bz and Eq^- In every case a 5-f unction initial distribution was used. In Figure 3 
Bz — lO^^T and Eqz ranges from 0.1 Vm~^ to 10 Vm~^. Not surprisingly, protons arc 
accelerated to higher energies as Eoz is increased. In contrast to the results obtained by 
Mori, Sakai, and Zhao (1998), the shape of the accelerated proton spectrum depends 
strongly on the model parameters: although in some cases the protons have a power 
law spectrum at low energy (1-100 keV for Eqz — 1 Vm~^), the power law index varies 
considerably. It should be noted that the frequencies / plotted in Figures 3 and 4 
represent the number of particles at the end of each simulation with energies in a fixed 
range of values of logio^^, i-e. 



where F is the true energy distribution. Thus, F{E) cx f{E)/E. At low energy, the 
distributions in Figure 3 correspond to F{E) oc E"^ where 7 < 2: Mori, Sakai, and 
Zhao (1998) obtained 7 ~ 2.0 - 2.2. 

Another striking feature of the distributions in Figure 3 is the formation of a bump- 
on-tail at high energy. This also occurs in the simulations of Petkaki and MacKinnon 
(1997), but not in those of Mori, Sakai and Zhao (1998). In general, models in which 
there is preferential acceleration of a small sub-population of protons to MeV energies 



f{E)d{log,oE) oc F{E)dE, 



(20) 
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and above are more efficient than models which predict monotonic decreasing spectra, 
since it is only at high energy that direct observational evidence for accelerated protons 
exists. The occurrence of bumps-on-tail in our simulations and those of Petkaki and 
MacKinnon, and their absence from those of Mori, Sakai and Zhao (1998), appears to 
be due to the choice of field configuration used in these studies. Mori and co-workers 
assumed a purely hyperbolic magnetic field in the {x, y) plane (i.e. 8B — 0) and a 
uniform E^. With finite S^, there is then a large parallel electric field component 
£■11 = EzB^/B"^ throughout the computational domain, and all the test particles are 
susceptible to strong acceleration. In our case, E\\ is significantly reduced outside 
the critical radius R = Rc, due to the presence of a perturbation to the hyperbolic 
magnetic field {SB ^ 0) and also a fall-off in E^ associated with the spatial profile 
of the reconnecting mode eigenfunction [cf. Equation (18)]. Consequently, particles 
initially lying inside R — Rc are subject to stronger acceleration than those initially 
lying outside. 
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Fig. 3. Final proton energy distributions for Bz — 10~^T and three 
values of Eq^. 

Final proton distributions are shown in Figure 4 for Eq^ — 1 Vm~^ and four different 
values of Bz. Results are not shown for Bz ~ since in this case the number of particles 
accelerated to energies of more than 1 keV was negligible. It can be seen that protons 
are accelerated to progressively higher energies as B^ is increased - in the case of 
Bz — lO^^T, particle energies of up to 40 MeV are observed. This appears to be due 
simply to improved particle confinement close to the X-point: the amplitude of the 
sinusoidal variation in Figure 1(a) varies as 1/Bz and the drift speed varies as l/-^^. 
For Bz < 10~^T bump-on-tail formation at high energy is again apparent. 
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Fig. 4. Final proton energy distributions for Eq^ — 1 Vm~^ and four 
values of B^. 

5 Discussion and Conclusions 

The addition of a longitudinal field component to a two-dimensional reconnecting 
magnetic configuration massively increases the efficiency of particle acceleration in such 
configurations. Physically, this is due to the fact that particles close to the magnetic 
X-point are strongly magnetized by the longitudinal field and arc not subject to strong 
grad-i?, curvature or E x B drifts. The effect of a small but finite longitudinal field Bz 
on particle acceleration is so dramatic that B^ = should be regarded as a singular case, 
unhkely to be representative of the conditions prevailing in the flaring corona. MeV 
protons can still be produced when B^ — 0, but the number of particles accelerated 
to such energies is only significant if the resistivity rj is assumed to be extremely large 
- much larger even than the anomalous values corresponding to ion acoustic or lower 
hybrid turbulence (Petkaki and MacKinnon, 1997). A key result of the present study 
is that the presence of a longitudinal guide field makes it possible for protons to reach 
MeV energies in a plasma with realistic values of rj. 

In the case of the simulation with B^ = 10~^T, Eqz — 1 Vm~^, 14% of the protons 
reach energies exceeding 1 MeV in Is. However, as discussed in Section 4.1, only a 
fraction of the initial particle population is modelled using the coUisionless equations 
of motion, and so the simulation results imply a super-MeV proton fraction of much less 
than 14%. As noted in Section 1, observations imply that up to 1-10% of the protons in 
the flaring corona are accelerated per second. While the test particle method adopted 
here yields valuable insights into the physics of proton acceleration in flares, a more 
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self-consistent approach may be required to meet the tight constraints on energetic 
proton fluxes imposed by recent observations. 

The sensitivity of the final distributions to (Figure 4) suggests that it might 
be possible to set constraints on the value of this parameter using 7-ray observations. 
However, it is not clear to what extent the simulation results depend on other features 
of the model, such as the assumption of an azimuthally-symmetric, centrally-peaked 
perturbation to the magnetic flux. CUEBIT could be used to study particle acceleration 
in X-point field configurations with a range of different normal mode perturbations. 
We also intend to study electron acceleration: the numerical scheme in Equations (2) 
and (3) allows sufficiently large timesteps (relative to the Larmor period) to make this 
feasible, although it will probably be essential to use a fully relativistic version of the 
code in this case. Collisions can also be added to the scheme in a straightforward way, 
enabling simulations to be carried out on longer timescales. 

In analytical studies of particle acceleration at reconnecting current sheets Litvi- 
nenko and Somov (1993) and Litvinenko (1996) formulated expressions for critical 
values of 5^, such that particles are efficiently accelerated, and for the typical energies 
which they can reach. Our results are not directly comparable to those of Litvinenko 
and Somov, since we have invoked a different field configuration. However, our par- 
ticle code could be applied to the particular current sheet geometry invoked by these 
authors and their analytical results compared with simulations. 
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